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Some models of diffusion-limited reaction processes in one dimension lend themselves to exact 
analysis. The known approaches yield exact expressions for a limited number of quantities of in- 
terest, such as the particle concentration, or the distribution of distances between nearest particles. 
However, a full characterization of a particle system is only provided by the infinite hierarchy of 
multiple-point density correlation functions. We derive an exact description of the full hierarchy of 
correlation functions for the diffusion-limited irreversible coalescence process A + A — > A. 
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Non-equilibrium kinetics of diffusion- limited reactions has been the subject of much recent interest ■ In contrast 
to equilibrium systems — which are best analyzed with standard thermodynamics — or reaction-limited processes — 
whose kinetics is well described by classical rate equations |^,^| — there is no general approach to non-equilibrium, 
diffusion-limited reactions. Some diffusion-limited processes in one dimension can be approached analytically, but 
the solution methods yield only a limited number of quantities of interest, such as the particle concentration, or 
the distribution of distances between nearest particles. In this letter, we derive an exact analytic description of the 
full hierarchy of multiple-point density correlation functions for the irreversible diffusion-limited coalescence process 
A + A — ► A. The hierarchy of correlation functions characterizes the system completely. 

Our model [9|-[ll|l is defined on the line — oo < x < oo. Particles A are represented by points which perform unbiased 
diffusion with a diffusion constant D. When two particles meet they merge into one particle which continues diffusing 
with the same diffusion constant D as the reacting particles. Since the reaction step is infinitely fast, the system 
models the diffusion-limited coalescence process A + A — > A. 

An exact treatment of the problem is possible through the method of Empty Intervals, known also as the method 
of Inter-Particle Distribution Functions (IPDF). The key parameter is E(x,y;t) — the probability that the interval 
[x, y] is empty (contains no particles) at time t. Particles just at the edge of an empty interval may diffuse into or 
out of the interval, affecting the probability E. With this observation in mind, one can write down a rate equation 
for the empty interval probability PJTc|]. For the model at hand, the rate equation is 

dE{x,y;t) d 2 d 2 

^i— = D{ dx^ + W ) iX ' m) - (1) 
The coalescence reaction imposes the boundary condition 

lim E(x,y;t) = 1, (2) 

y\x or xly 

and in addition, as long as there are any particles left in the system, 

lim E(x, y;t) = Q. (3) 



From E(x,y;t) one can derive useful parameters, such as the concentration of particles; 

dy 

or the probability that given a particle at x the next nearest particle is at y (the IPDF) 



p(x;t) = -— E(x,y;t)\ y=x , (4) 



d 2 

P{x,y,t) = p(x,y;t) -g^^ E(x,y;t) . (5) 



While such level of description affords us invaluable physical insights into many interesting applications 1 1 1 , it 



does not characterize the system completely. For that we need the multiple-point density correlation functions 
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p n (x±, X2, ■ ■ ■ , x n ; t), i.e., the joint probability density to find n particles at xi, X2, ■ • • , x n at time t. For n = 1, pi(x, t) 
is identical with the particle concentration of Eq. However, a complete characterization of the system requires 
knowledge of the full hierarchy of correlation functions {p n }n°=i- 

The multiple-point correlation functions may be obtained from a generalization of the method of Empty Inter- 
vals pc[| . Let E n (x\,yi,X2,y2i ■ ■ ■ )X n ,y n ;t) be the joint probability that the intervals (i = 1,2, ... ,n) are 
empty at time t. The intervals are non-overlapping, and ordered: x\ < y\ < ■ ■ ■ < x n < y n . Then, the rt-point 
correlation function is given by 

gn 

p n (xi, ...,x n ;t) = {-l) n -^——g—E n (x 1 ,y 1 , . . . , x n ,y n ; t)\ m = Xl ,..., Vn = Xn . (6) 
Doering [fLOl has shown that for irreversible coalescence the E n satisfy the partial differential equation: 



d d 2 d 2 d 2 d 2 

-E n {x 1 ,y 1 ,...,x n ,y n -,t) = D(^ + w ^--- + —J—)E n , (7) 



with the boundary conditions 



and 



lim E n (xuVi, ■ ■ - ,x ni y n ;t) = E n ^ i (x 1 ,y l , ... , fa, fa, . . . , x n , y n ; t) , (8) 

XilVi or yil'Xi 



lim E n (xi,yx, . . . , x n ,y n ;t) = E n -i(x\,yi, . . . , fa, fa+i,. ..,x n , y n ; t) . (9) 

yi]x i+ i or x i + ilyi 



For convenience, we use the notation that crossed out arguments (e.g. $i) have been removed. The E n are tied 
together in an hierarchical fashion through the boundary conditions (g) and (|)|): one must know E n ^\ in order to 
compute E n . 

We now provide with a full solution of the E n . At the level of n — 2, the solution is |l2[ : 

E 2 (xi,yi,X2,y2;t) = 

E(x 1 ,y 1 ;t)E(x 2 ,y2;t) - E(xi, x 2 ; t)E{y u y 2 ; t) + E(x 1 ,y 2 ;t)E(y 1 ,x 2 ;t) . (10) 

Clearly, each term on the r.h.s. of Eq. ( |lO| ) satisfies the equation ([?]), as well as the boundary conditions (||) and (|^) 
— since E(x, x; t) = 1 (Eq. ||). Similarly, E n can also be expressed in terms of products of E's for single intervals: 

(2n-l)!! 

i yi ) • ■ ■ j x n , yn ;t) = a p E(zi < p,Z2, p ;t)E(z3 ! p,Z4,p;t) ■ ■ ■ E(z2 t) . (11) 

p=l 

Here Z\^ p , Z2 tP , . ■ ■ , Z2 n ,p symbolize an ordered permutation, p, of the variables x\, y%, . . . , x n , y n , such that 

%l,p < Z2,p, Zz iP < Z4 iP , . . . , Z2n-l,p < Z2n,p, and < Z3 iP < 25 iP • • • < Z2n-l,p ■ (12) 

There are exactly (2n — 1)!! = 1 • 3 (2n — 1) such permutations. a p is +1 for even permutations (permutations 

that require an even number of exchanges between pairs of variables), or —1 for odd permutations. 

Eq. (0) can be proved by induction. For n = 1 it reduces to Ei(x,y;t) = E(x,y;t), and for n = 2 it reduces 
to Eq. flic]), as required. We need only show that if (10) is true for n — 1 (n > 3), then it is also valid for n. It is 
easy to see that because each of the E's satisfies Eq. (|l]), the proposed E n satisfies Eq. (^). Now test the boundary 
conditions: Suppose that Xi = yi. The permutations in the r.h.s of (|ll| ) are divided into two groups: (a) those which 
leave (xi,yi) as an argument of one of the E's in the product, and (b) those which separate Xi and yi into different 
E's in the product. The permutations in group (a) add up to E n _i(xi, y%, . . . , $i, fa, ... , x n , y n ; t), as required by 
the boundary condition ( ®). This is because E(xi,yi\) = 1 (from Eq. ||), and because of the induction assumption 
regarding the validity of ( |ll| ) for n — 1. Notice that the parity of the permutations in group (a) is the same as if Xi 
and yi were removed. On the other hand, the permutations in group (b) add up to zero, for the following reason. 
Suppose that in some permutation Xi and yi are paired with other variables z\ and Z2, (x{,z\) and (j/j,^)- Then, 
there exists a similar permutation of the variables x\, y\, ■ ■ ■ , x n , y n where individual pairs remain in the same order, 
but now the pairings of Xi and yi are exchanged; (xi, Z2) and (yi, z\). When Xi — yi, the product of the E's in these 
two permutations is identical. But the parity of the two permutations is opposite, and so they add up to zero. The 
same is true if the pairings of Xi and yi are (zi,Xi) and (z2,yi), or (z\,Xi) and (yi,Z2). [Notice that the pairings 
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(xi,z\) and [z2,Vi) cannot occur, because of the required ordering, Eq. (|T^).] In sum, the boundary condition (||) is 
satisfied. The proof of @ follows a similar line of reasoning. 

The ordered permutations of the endpoints of n intervals (2n variables) may be constructed recursively in the 
following way. The order constraint of Eq. (|l^ ) requires that zi_ p — x\, for all permutations p. Set then Z\ = x\, and 
zi equal to one of the other yi, . . . ,x n ,y n variables. Then, arrange the remaining 2n — 2 variables in all their ordered 
permutations. Finally, repeat this procedure, selecting sequentially z-i — yi, X2, V2, ■ ■ ■ , x n , y n . Thus, the number 
of permutations for the endpoints of n intervals, N(n), satisfies the recursion relation N(n) = (2n — l)N(n — 1), 
and therefore N(n) = (2n — 1)!! — since iV(l) = 1. The recursive construction allows us also to express E n more 
compactly, in terms of E n -\. 

n 

t) = +J2 E{x 1 ,y j -t)E n _ 1 {?t ll y l , . . . , Xj , fa, . . . , x n ,y n ; t) 

3 = 1 

Tl 

-^2E(x 1 ,x j :t)E n - 1 (tf 1 ,y 1 , ^j,yj, . . .,x n ,y n ;t) . (13) 
i=2 



Until now we have ignored the issue of initial conditions . For the solution ( |11| ) to work, it is required that the same 
relation be satisfied at time t = as well! This seems at first sight a formidable restriction, but fortunately some most 
important situations are unaffected by it. If the particles are initially randomly distributed, independently from each 
other, at a homogeneous concentration po, then E(x,y;0) = e~ p °( y ~ x ' > and 

E n (x 1 ,y 1 ,...,x n ,y n ;0) = e-r°^-^ + -+^-^ . (14) 



One can easily check that Eq. ( |ll| ) is satisfied. In this case, Eq. (|6|) yields the obvious relation: p(x%, X2, ■ ■ ■ , x n ; 0) = Po • 
This uncorrclatcd random initial distribution gets quickly correlated with time. 

Another interesting situation is the state of the system in the long-time asymptotic limit. It can be shown that 
the system arrives at a universal asymptotic state, independent of the initial distribution of particles. (This excludes 
some exotic situations, such as fractal initial distributions.) The long-time asymptotic solution of Eq. (|I]), with the 
boundary conditions (||) and (||) is 

E(x,y;t) = erk(2^JL) . (15) 
Using Eq. (Q), the corresponding long-time asymptotic concentration is: 

PasympM = vibr (16) 



From Eqs. ( pi] ) and (^|) we get the two-point correlation function 



P2( ^ 2;t) = 1 - e-^ + V^ 2 erfc(C) , (17) 

/^asymp. 

where we used the notation £ = (#2 — xi)/\8Dt. For the three-point correlation function, we get 

p 3 (x u x 2 ,x 3 ;t) = 1 _ e _ 2c2i _ e _ 2& _ e _ 2?|i + 2e -g 1 -e 32 -& 

Pasymp. 

+^61 (e^ 1 - e _,5 ' 2_ ^i)erfc(^ 2 i) 
+Vn £,32(e~ (?i2 - e-^i-«!i) er fc(£ 32 ) 

+^&i(e~& - e"^ _ ^ 2 )erfc(^ 31 ) , (18) 

where now ^ = (xi — Xj)j\J 8Dt (notice that £31 = £32 + £21 is not an independent variable). 

In Fig. [I], we show the two-point correlation function in the long-time asymptotic limit (Eq. [l?]). We see that the 
two points become uncorrelated as the distance between them increases, but that there is an effective strong repulsive 
interaction (due to the calescence reaction) between nearby particles. Interestingly, the two-point correlation is a 
monotonous function of the distance. A simple convolution of the distances between nearest particles predicts an 
oscillating tail p3[. 
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The three-point correlation function, Eq. (|lq), is a bit harder to illustrate. Instead of a full description, in Fig. g 
we compare p$(x\, X2, X3) to P2(%i,%2)p2(x2,%3)/p(%2) (in the spirit of the truncation ansatz that might be used in a 
Kirkwood approximation, say), along the line x% — a; 2 = xi — X\. Again, we see that as the distance between the three 
particles increases they become rapidly uncorrelated, but that the approximation ansatz fails for short distances, due 
to reactions. 

In summary, we have obtained the complete hierarchy of 71-point density correlation functions for the diffusion- 
limited irreversible coalescence process, A + A — ► A, in one dimension. The recursive form of Eq. ( plf ) should be useful 
to prove theorems. Moreover, with today's available software for symbolic manipulation, it is relatively easy to churn 
out explicit expressions for a finite number of points. In fact, we have used "Mathematica" to compute the examples 
of 2- and 3-point correlation functions. 

The more general problems of reversible coalescence — when the back reaction A — > A + A is allowed — and 
coalescence with different kinds of particle input, can also be handled, in principle, by the empty interval approach. 
That is, in both cases (even both at the same time) there is a closed hierarchy of linear partial differential equations, 
coupled through their boundary conditions, similar to Eqs. (0), (||) and @ JlOj . But the linear operators involved 
break the simple x-y symmetry of Eq. (0), and a complete solution remains an open challenge. 

ACKNOWLEDGMENTS 

I thank Larry Glasser for useful discussions, and Charlie Doering for sharing with me his elegant solution for the 
two-point correlation function. 



[1] N. G. van Kampen, Stochastic Processes in Physics and Chemistry (North-Holland, Amsterdam, 1981). 
[2] H. Haken, Synergetics ( Springer- Verlag, Berlin, 1978). 

[3] G. Nicolis and I. Prigogine, Self-Organization in Non-Equilibrium Systems (Wiley, New York, 1980). 
[4] T. M. Liggett, Interacting Particle Systems (Springer- Verlag, New York, 1985). 

[5] K. Kang and S. Redner, Phys. Rev. A 32, 435 (1985); V. Kuzovkov and E. Kotomin, Rep. Prog. Phys. 51, 1479 (1988). 
[6] J. Stat. Phys. Vol. 65, nos. 5/6 (1991); this issue contains the proceedings of a conference on Models of Non-Classical 

Reaction Rates, which was held at NIH (March 25-27, 1991) in honor of the 60th birthday of G. H. Weiss. 
[7] K. J. Laidler, Chemical Kinetics (McGraw-Hill, New York, 1965). 

[8] S. W. Benson, The Foundations of Chemical Kinetics (McGraw-Hill, New York, 1960). 
[9] D. ben-Avraham, M. A. Burschka, and C. R. Doering, J. Stat. Phys. 60, 695 (1990). 
[10] C. R. Doering, Physica A 188, 386 (1992). 

[11] For reviews, see D. ben-Avraham, Mod. Phys. Lett. B 9, 895-919 (1995); D. ben-Avraham, in Nonequilibrium Statistical 

Mechanics in One Dimension, V Privman, ed., pp 29-50 (Cambridge University Press, 1997). 
[12] M. A. Burschka, C. R. Doering, and W. Horsthemke, unpublished. 
[13] P. A. Alemany and D. ben-Avraham, Phys. Lett. A 206, 18 (1995). 

FIG. 1. Two-point correlation function for the coalescence process in the long-time asymptotic limit. Shown is P2(£)//°Lymp. 
vs. f. 

FIG. 2. Three-point correlation function for the coalescence process in the long-time asymptotic limit. Plotted is the relative 
error made by the Kirkwood approximation, (pxirkwood — pz)l 1 pz, for the line £21 = £32 = £■ The Kirkwood approximation in 
this case is pxirkwood = P2(£) 2 /p asymp. • 
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